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Nowadays the sector decomposition technique, which can isolate divergences from parametric 
representations of integrals, becomes a quite useful tool for numerical evaluations of the Feynman 
loop integrals. It is used to verify the analytical results of multi-loop integrals in the Euclidean 
region, or in some cases practically used in the physical region by combining with other methods 
handling the threshold. In an intermediate stage of the sector decomposition for the multi-loop 
integrals, one often has to handle enormously large expressions containing tons of terms. The 
symbolic manipulation system FORM is originally designed to treat such huge expressions and has 
a strong advantage for it. In this talk, the implementation of the sector decomposition algorithm 
on FORM is discussed. A number of concrete examples including cases of multi-loop diagrams 
are also shown. 
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1. Introduction 

The Large Hadron Collider (LHC) has started its operation, and International Linear Collider 
(ILC) is planned as a next generation collider. In order to obtain as much information as possible 
from such experiments, it is undoubtedly important to know the precise predictions from the theory 
at high energies. This fact pushes forward recent study on the computation of one- and higher loop 
corrections in the perturbation theory. 

To evaluate the loop integrals, one may try to find formulae of them analytically. However, if 
they have many loop momenta, many external legs, and many kinematic parameters, it becomes an 
extremely difficult problem. Another approach to the evaluation of the loop integrals is to perform 
loop integrals by using a numerical integration method. In principle, it is possible to construct 
a general recipe for computing loop integrals numerically, which is applicable to a wide enough 
class of the loop integrals. Even if one has an analytical formula for some integral, it is good idea 
to have another method for computing any integrals numerically, for checking purpose. In practice, 
to perform loop integrals numerically, one has to handle several singularities: UV, IR divergences, 
and singularities from physical thresholds. 

In this talk, the implementation of the sector decomposition algorithm [[T], ^ 1 which can isolate 
divergences from parametric representation (e.g., Feynman parameters) of integrals, on a symbolic 
manipulation system FORM is discussed. Then, an application of the numerical extrapolation 
method on handling physical thresholds [[7|] after isolating all IR divergences is also discussed. 

2. The sector decomposition for IR divergent loop integrals 

Here we briefly sketch how the sector decomposition algorithm is applied for IR divergent 
loop integrals, in a easy example. We consider the one-loop scalar massless on-shell box diagram 
in D = 4 — 2e dimensions. By introducing Feynman parameters and performing the momentum 
integration, we arrive to 



where s = {p\+ P2) 2 , t = (p2 + P3) 2 , and pf = (i = 1, . . . ,4), and we have used abbreviations, 
P\2 = Pi + P2, P123 = Pi + P2 + P3 and X1234 = xi +X2 + X3 + X4. This integral has overwrap sin- 
gularities on edges of the integration region: the denominator of the integrand vanishes when two 
variables simultaneously go to zero, e.g., x\ ,X2 — > 0. The sector decomposition algorithm can disen- 
tangle such overwrap singularities, by splitting integral domain and making variable replacements, 
iteratively. One can obtain 
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<ee [O] for a review. Currently, available public codes of the sector decomposition are M, ra 
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-s - tx 2 x 3 - i0) 2+e 



+ (s<->t). (2.2) 



Note that, in this form, there are no overwrap singularities, and all singularities are factorised as 
powers of monomials. This can be expanded with respect to £ as: 

7 4 0m /r(2 + £) = p + ^ + C + ^(£), (2.3) 



with 



r l i 

C 2 =4/ dx- —2 + (j«->f), (2.4a) 
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Ci = -4 / dx-± — 1 - 2 / d 2 x- — 7 + (s++t), (2.4b) 
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The coefficients of Laurent expansion C 2 , C\ and Co are now expressed in terms of multidimen- 
sional integrals, and in the Euclidean region (s < and t < 0), they can be easily computed. 

3. Implementation of the sector decomposition on FORM 

In an intermediate stage of the sector decomposition for more complicated loop integrals, one 
needs to handle a lot of terms and very large expression. For example, the number of the generated 
sub-sectors for the scalar massless on-shell triple box diagram is ^"(10000), and then many more 
terms are produced by the £-expansion. Most of popular computer algebra systems try to keep 
the expressions in the physical memory. When the expressions become large, the disk memory is 
used via the virtual memory manager of the computer, and then a computer algebra system can 
be extremely slowed down or does not work. Therefore, in practice, there is a limit of the size of 
expressions in such systems 2 . FORM is designed to treat such huge expressions, which are larger 
than the available physical memory, by using the disk with less penalties in the performance, and 
therefore has strong advantage. 

Our program uses FORM for symbolic manipulations. From the user input file describing 
the integral, FORM produces Fortran code for the integrand, after the sector decomposition. The 
program has been checked by various multi-loop integrals, which include planar and non-planar 
massless on-shell double boxes up to €?(e°), massless three-loop propagators up to i^(£ 2 ), etc. As 
an example, the numerical result of the planar massless on-shell double box is shown in Table [|. 

4. Handling physical thresholds by the numerical extrapolation 



Consider the one-loop box diagram Eq. (2.3) again. If either s or t is not negative and the 



denominator (or the argument of the logarithm) of the integrand is not positive definite, it can be- 
come to zero at some points in the integration region. For such singularities arising from physical 
2 To avoid this limit, FIESTA ^ uses a data base manager for storing the expressions on the disk efficiently. 
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Result and Error 


Elapsed Time 


c 4 


0.2000000066 1904E+01 ± 0. 124855 17866572E-05 


<0.1s 


c 3 


-0.6000000 1069 17 1E+01 ± 0.540967867 12553E-05 


<0.1s 


c 2 


-0.49167420522766E+01 ± 0.17982248269768E-04 


0.3s 


Ci 


0.11494738132380E+02 ± 0.85267241358588E-04 


9.5s 


Q, 


0.13801183392483E+02 ± 0.299801 15734943E-03 


134.6s 



Table 1: The numerical result of the planar massless on-shell double box for s = — 1 and t — — 1 (CPU: 
Xeon 5160 3GHz). C,- are the coefficients of 1 /e' (overall T(3 +2e) is excluded). With this parameter set, 
the integrand is relatively smooth (no strong peaks), and we used DCUHRE as the multidimensional 
integrator. The goal of relative tolerance of the numerical integration is set as E Te i = 10~ 6 for each integral. 
The boldfaced digits indicate they agree with the analytical values. 



thresholds, the contour deformation of Feynman parameters is used in Actually the com- 
bination of the sector decomposition and the contour deformation is applied in some practical 
calculations [10, 11]. 

Here we consider another possibility, the numerical extrapolation method [^. In this method, 
we put i8 with a small but finite 8 instead of iO in the denominator. Then we can compute the 
integral 1(8) for a given 8. Calculating the sequence {I (8k)} for 8k = 5ir* _1 (k = 1,2,3,..., 
< r < 1) and extrapolating them by an adequate method, we can finally obtain the result of the 
integral in the limit of 8 — ► 0. 

The numerical result of the massless one-loop on-shell box and the massless one-loop box 
with two adjacent off-shell legs are shown in Table || and ||, respectively. 



5. Summary 

The implementation of the sector decomposition algorithm on a symbolic manipulation sys- 
tem FORM is discussed. Thanks to FORM'S advantage for handling very large expressions, the 
limitation due to the amount of available physical memory is greatly relaxed. After isolating all IR 
divergences by the sector decomposition, the numerical extrapolation is used for handling threshold 
singularities. Several numerical results are also shown. 
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Result and Error 


Elapsed Time 


c 2 


real 
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0.1s 


Co 


real 
imaginary 
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0.8s 
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Result and Error 


Elapsed Time 


c 2 


real 
imaginary 
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real 
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